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Abstract. We propose a two-dimensional phase-field-crystal model for the (2x1)- 
(1x1) phase transitions of Si(OOl) and Ge(OOl) surfaces. The dimerization in the 
2x1 phase is described with a phase- field-crystal variable which is determined by 
solving an evolution equation derived from the free energy. Simulated periodic arrays 
of dimerization variable is consistent with scanning-tunnclling-microscopy images of 
the two dimerized surfaces. Calculated temperature dependence of the dimerization 
parameter indicates that normal dimers and broken ones coexist between the 
temperatures describing the charactristic temperature width of the phase-transition, 
Tl and Th, and a first-order phase transition takes place at a temperature between 
them. The dimerization over the whole temperature is determined. These results 
are in agreement with experiment. This phase- field-crystal approach is applicable to 
phase-transitions of other reconstructed surface phases, especially semiconductor nxl 
reconstructed surface phases. 
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1. Introduction 

Semiconductor surfaces are of huge importance especially in the modern era of 
nanoscience and nanotechnology. Usually, a bulk terminated surface (1 x 1) is unstable, 
will undergo a surface reconstruction and become a stable reconstructed surface phase 
(m X n) P, El [3]. Most of reconstructed surfaces transit to some 1x1 structures 
at elevated temperatures. The Si(OOl) and Ge(OOl) reconstructed surfaces have been 
extensively studied because they are closely relevant to the modern computer technology 
[H [5]. For both of them, the reconstruction is realized through the forming of 
regular arrays of dimers in the top layer (2 x 1 or dimerized phase), which has been 
confirmed by scanning-tunnelling-microscopy (STM) experiment [HIUlS]. It was shown 
experimentally [6l [7] that when heated to certain temperatures, the 2x1 reconstructed 
surfaces will transit to 1 x 1 structures. There were some first-principles calculations for 
the local atomic configurations of the 2x1 surfaces [8], but it is controversial even for 
the essence of the phase transitions [6l[7j. Considering that a structural phase means an 
averaging of local atomic structures over large enough scales, the reconstructed surface 
phases need further clarification and the essential physics of the phase transitions are 
still unknown. A theory that can elucidate the issues is highly desirable. 

Phase-field method is a reliable approach to modeling and simulating structural 
phases and dynamical phase transitions P [lO]. It has been applied to various fields 
such as spiral surface growth [11], dendritic growth [12], alloy sohdification [13], crystal 
nucleation [Tl], step- flow growth [15], [IE], epitaxial island growth [13, [181 [I9] , and surface 
phase transition dynamics [20]. Recently, a phase-field-crystal approach was proposed 
to model the internal spatial structures of a given phase [21], and periodic lattices was 
obtained by solving evolution equations derived from the free energies. This powerful 
approach has been successfully used for natural modeling of elastic interactions [22] and 
binary alloy solidification [23] . 

In this paper we propose a phase-field-crystal model for the Si(OOl) and Ge(OOl) 
dimerized surface phases and their phase transitions. We use a two-dimensional (2D) 
phase-field-crystal variable to describe the dimerization of atoms in the top layer, and 
determine the variable by solving an evolution equation derived from the free energy. 
Our simulated morphology of periodic arrays of the dimerizarion variable is in good 
agreement with large-scale STM images of the Si(OOl) and Ge(OOl) surfaces in high- 
quality samples [21 [3l HI [5]. Furthermore, we derive the temperature dependence of 
the dimerization parameter, and show that normal dimers and broken ones coexist 
between two characteristic temperatures, Tl and Th, and a first-order phase transition 
takes place at a phase-transition temperature Tc in between Tl and Th- These are in 
agreement with experiment [6l [7] . This phase- field-crystal approach can be applied to 
other reconstructed surface phases and their phase-transition dynamics. 

The remaining part of this paper is organized as follows. In next section we shall 
present our phase- field-crystal modeling for the 2x1 phase and its phase transition 
dynamics during transiting to the 1x1 phase. In section III we shall present our main 
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simulated results. In section IV we shall apply the model and the simulated results to 
the (2xl)-(lxl) phase transitions of Si(OOl) and Ge(OOl) dimerized surfaces. Finally, 
we shall give our conclusion in section V. 

2. Phase-field-crystal modeling 

There are two kinds of the (001) surfaces, type-S^i and type-SB, for both Si and Ge. 
They appear alternately in the vertical direction, but do not exist in the same layer in 
high-quality samples, although they can be changed into each other by a 90° rotation 
[H [5]. Because the dimerization takes place only in the top layer, it is reasonable to 
describe both of the Si(OOl) and Ge(OOl) dimerized surface phases by a 2D model. 
We neglect the buckling of dimers without losing main physics because we are mainly 
interested in their high-temperature phase transitions to corresponding 1x1 phases. 
We show the atomic configurations and corresponding phase-field-crystal modeling in 
Fig. 1. We use r = for a 2D point, where x describes the coordinate in the 

dimerization direction or the horizontal direction and y in the other direction. The 
key order parameter is the dimerization-induced change of the bond length, Ab, in the 
horizontal or x direction because there is no dimerization in the other direction P, [7j. 
Phase-field-crystal variable (j){x, y) is periodic in the x direction and uniform in the other 
direction. Along the x direction, y) reaches the maximum, 0max, at the center of 
the dimer and the minimum, — 0max, at the middle point of the two nearest dimers with 
the same y. This model is consistent with STM images of regular arrays of dimer chains 
in high-quality Si(OOl) and Ge(OOl) dimerized surfaces [H El El [7j . The average value 
of y) over r, 0, is set to zero for the dimerized surfaces. The variable y) for 
the dimerized surfaces is normalized so as to make the integration of \(f){x,y)\ over the 
2D unit cell equal 2A6 because there are two top-layer atoms in the 2D unit cell. The 
variable y) is identical to zero for the undimerized surfaces. These are enough for 
the following phase-field simulation, but, with all these in our mind, we can imagine 
that the morphology of the phase-field variable y) can be described approximately 
by a simple function cos(7ra;/c), where c is the lattice constant of the undimerized 
2D unit cell and the center of the dimer is taken as the zero point of x. 
The free energy of </)(x, y) can be written as 

T = \j^ dr{\ct>[A + {ql + - + «6^}, (1) 

where the A is defined hy A = a(T — Tl), and the parameters A, u^, uq, and a are 
positive constants, go is defined as 27t/cq and cq has the meaning of lattice constant if 
periodic solution is obtained. The operator G = (gp + is constructed by fitting to 
the dimerization structure |2T]. The M4 and uq terms, as usual, are used to describe the 
first-order phase transitions in the Si(OOl) and Ge(OOl) surfaces. The bilinear term is 
sued to describe the effect of temperature and dimerization structure. The parameter uq 
can be set to 1 by redefining the phase-field-crystal variables without losing any physics. 
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Figure 1. The phase-field-crystal modeling of the Si(OOl) and Ge(OOl) dimerized 
phases, (a) shows atomic models of the dimerized (001) surfaces and (b) that of 
bulk-terminated (001) ones. The dimerization happens only in the horizontal (or x) 
direction. The gray circle is the atom in the top layer, the bigger black dot in the 
second layer, and the smaller black dot in the third layer. The curve (c) is schematic 
phase-field-crystal (j) description of the dimerization-induced atomic density change 
with respect to that without dimerization (d). 



The evolution equation of is a time-dependent Cahn-Hilliard (CH) equation [21] 

where 77 is a Gaussian random variable whose average value is set to zero. The first 
term on the right side describes the diffusion of the phase-field variable. Using the 
dimensionless phase- field variable ip, dimensionless coordinates x, and dimensionless 
time r defined by [21] 



X = (x, y) = qor, ^ = ^^4>, r = XTq^t, (3) 
we express the free energy (1) and the CH equation as 



and 



^^^„/^.f{||e + (l + -)>-^V.' + -V.»} (4) 

^ = V2{|e + (l + ^n^-«^» + ^^} + C (5) 

Now we have only three independent parameters: jFg, e, and u. has the dimension of 
energy and the other two are dimensionless. They can be expressed in terms of original 
parameters: Tq = Xq^j^JuQ, £ = a{T — TLj/q^ = aiT/Ti — 1), and u = u^/ {q^y/uo). 
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The new random variable ( whose average value is also zero is subject to the two-point 
correlation function 

{({x, t)C{x', t')) = DV^6{x - x')6{t - t') (6) 

with D = ksT/J^Q. The average value of ip{x,y), ip, is conserved by Eq. (5) and 
therefore ip can be taken as an independent parameter. 

3. Main simulated results and analysis 

We solve the dimensionless evolution equation Eq. (5) by difference method. The 
Laplace operator and the second-order differential operator d'^/dx'^ are discretized 
by the central difference formula as usual, and the time coordinate is discrerized using the 
first-order finite differential approximation. The periodic boundary condition is adopted 
for all the simulations. The initial condition is subject to Gaussian random fluctuations. 
We find that resultant equilibrium patterns are independent of the random variable (. 
We keep ip = because the high-temperature 1x1 phase has ^ = 0. The parameter 
JFq is not directly relevant in our solving Eq. (5). We have tried various values for the 
system sizes and the parameters u and e. We choose approximately 40 periods or 80 
surface lattice constants when doing the main simulations, but confirm our results using 
larger systems. We use 6x = vr/lG and St = 10~^ for the space and time increments. 
The morphology presented in the following is only a part of the whole system. 

Presented in Fig. 2 is the periodic morphology of phase-field-crystal variable ip{x, y) 
with the parameters: u = 0.1, e = —0.005, and ip = 0. The period is equivalent to 27r 
within an error of 10~^. The center of the white stripe corresponds the maximum 
ipma.x and that of the black stripe the minimum ipmin- This phase- field-crystal simulated 
dimerization pattern is in agreement with experimental STM images of high-quality 
Si(OOl) and Ge(OOl) dimerized surfaces [21 [3], HI [5] . This stripe pattern is obtained for 
the parameter region: 0.00221 > e > —0.01. 




Figure 2. Simulated striped morphology of phasc-field-crystal variable ijj (upper part) 
and its side view curve (lower part). We use the parameters Sx = 7r/16, St = 10^^, 
u = 0.1, e = -0.005, and ^ = 0. 
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Figure 3. The calculated parameter i? as a function of e = a{T — Tl), with the 
phase-transition point Be indicated by dash line. The e-increasing part is shown 
by the solid line with circles, and the e-decrcasing part by the solid line with 
squares. The phase coexistence region is defined by < e < eh- The parameter 
set {eltSc^h), corresponding to {Tl,Tc,Th), is (0.0,0.169,0.221) for m=1.0, and 
(0.0,0.00169,0.00221) for u^O.l. 



We use B = ipmsix — V'min OUT Order parameter and present its e dependence for 
u=1.0 and u=0.1 in Fig. 3. It is clear that the phase transition is of first order. The 
B, as a function of e, in the case of n = 0.1 is ten times smaller than that in the case 
of M = 1.0. The e width of the phase coexistence region, eh — £l, is proportional to 
M^. Because e is proportional to T — Ti, the e dependence implies the temperature 
dependence from zero temperature to Tl, and finally beyond Th- It is proved using a 
series of simulated results that Tc is also proportional to u"^. Therefore, the parameter 
u describes the temperature width of phase coexistence region. Systematical analysis of 
simulated results shows that ip can be quantitatively described by the function 

^ = ^sin(— x), (7) 
^ Qo 

especially in the neighborhood of the maximum and the minimum. Here B and q are 
determined by the simulated results. \q/qo — 1| can be very small as long as u is small 
enough. Actually, q is equivalent to go within a tiny error less than 10"^ as long as u is 
smaller than 0.1. 

On the other hand, B and q can be analytically determined by minimizing the free 
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energy (4) in terms of a variational expression ip = C sm{px) similar to Eq. (7). In this 
way, we derive p = 1 and 



1 



C = -^y3n+ ^9n2 -40a(T/TL - 1). (8) 

This means q = Qq. The expression (8) is reasonable only when T is not larger than T^, 
which implies that u can be expressed as 
2 



u = -^10a{TH-TL)/TL. (9) 

We can obtain an analytical expression of the order parameter B = 2C for n < 0.1 by 
requiring that C given by Eq. (8) is equivalent to B/2 obtained by numerically solving 
Eq. (5). The phase-transition temperature is given by 

Tc = + \tl. (10) 
Using the relation </> = qQ-ip (due to = 1) and Eq. (9), we obtain B = Ab/c and 

Ab = 2c 



for T < Th, and B = Ab = for T > Th- This expression determines the temperature 
dependence of the dimerization-induced change of the bond length, Ab. 



4. Applied to the (2xl)-(lxl) phase transitions of Si(OOl) and Ge(OOl) 

For the Si(OOl) and Ge(OOl) dimerized surfaces, it is relatively easy to measure Tl, Th, 
Co (zero temperature), and AbuT (room temperature, 300 K). It should be noticed that 
q = 2-71 /c, where c = 2c (cq = 2co). We use these parameters as input. The parameter a 
can be determined by substituting 300 K and A6rt for T and Ab in Eq. (11), and then 
calculate u and in terms of Eqs. (9) and (10). We calculate U4 using the definition 
M4 = uq^ (due to Uq = 1). The input and calculated results are summarized in Table 
I. In Fig. 4 we present the dimerization parameter Ab as a function of temperature for 
both the Si(OOl) and the Ge(OOl) surfaces. It is clear that Ab is still finite at T^, and 
jumps to zero once T is larger than Th- 

It should be pointed out that our phase-field-crystal equilibrium patterns, as shown 
in Fig. 2, consist of regular arrays of infinitely long dimer chains. These are in agreement 
with large-scale experimental STM images of the parallel perfect dimer chains in the 
high-quality Si(OOl) and Ge(OOl) surfaces [HIS]. Generally speaking, such a dimer chain 
can be broken into several segments, but these segments are still completely in the same 
line [H El El [31 El [7j and can be considered to be an infinitely long dimer chain in the 
sense of averaging along the line. Actually, there is either type-S^ or type-S'^ dimerized 
phases in high-quality (001) surfaces [HE]- Therefore, our phase-field-crystal theory can 
describe well the experimental regular arrays of dimer chains. 

In addition, our simulated results are in agreement with experimental observation 
that the dimers still exist above Tc [6], as shown in Fig. 4. Some of the dimers begin 
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Table 1. The input parameters and the calculated results for the dimerized Si(OOl) 
and Ge(OOl) surfaces. 









Co [A) 






1223 


1473 


3.84 


1.6 




a 


u 








5.8x10"^ 


0.022 


0.06 


1410 


Ge(OOl) 


Tl{K) 


Th{K) 


Co (A) 


A6rt (A) 




950 


1130 


3.99 


1.6 




a 


u 




Tc{K) 




5.4x10"^ 


0.020 


0.05 


1085 



1.8 




^0 300 600 900 1200 1500 

T(K) 



Figure 4. The temperature dependences of dimerization parameter Ab (A) for the 
Si(OOl) (square) and Ge(OOl) (circle) phases, with Tc indicated by dash line. Th is 
1473 K for the Si(OOl) and 1130 K for the Ge(OOl). Ab is finite at and jumps to 
zero once T > Th- 

to break at Tl and all of them finally disappear at Th- Between and T^, there is a 
phase coexistence of the normal dimers and the broken ones, but there is no coexistence 
of the type-S^ and type-Ss discussed in earlier references [25]. This is the clear sign 
of the first-order phase transition and thus there must be some dimers above the phase 
transition temperature Tc- Naturally, Th — Tc varies in different samples, and can be 
so small that Th is equivalent to Tc within measurement error, which leads to some 
conclusions that there is no dimers above Tc [7]. Essentially, a structural phase can 
be judged only when it has a large enough spatial size. Therefore, our modeling and 
simulated results are not only reasonable but also in agreement with experiment [6l [7] . 

5. Conclusion 

In summary, we have proposed a two-dimensional phase-field-crystal model for 
the Si(OOl) and Ge(OOl) dimerized surface phases and their phase transitions to 
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corresponding 1x1 phases at elevated temperatures. We use a phase-field- crystal 
variable to describe the dimerization of atoms in the top layer, and determine it 
by solving the evolution equation derived from the free energy. The simulated 
periodic arrays of dimer chains are consistent with STM images of the Si(OOl) and 
Ge(OOl) dimerized surfaces. The calculated temperature dependence of the dimerization 
parameter shows that normal dimers and broken ones coexist between and Tjj and 
the first-order structural phase transition takes place at in between. These results are 
in agreement with experiment. This phase- field-crystal approach can be directly applied 
to phase transitions of semiconductor n x 1 reconstructed surface phases, and should be 
suitable to other semiconductor reconstructed surfaces and their phase transitions. 
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